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Abstract. We present a quantitative study of the morphology of 41 Extremely Red 

Objects (EROs). The analysis is based on deep otical and near-infrared images from the 

Hubble Space Telescope public archive, and performed by fitting to each galaxy image a 

PSF-convolved bi-dimensional model brightness distribution. Relying both on the visual 

inspection of the data and on the results of the fitting procedure, we are able to determine 

q , the fraction of irregular and/or interacting EROs, and to identify those that more closely 

^ ■ resemble local ellipticals. To the former class, whose members arc probably high redshift 
C/2 . 

■ dusty starburst, belongs about 15% of the whole sample, whereas the elliptical-like ob- 
, ^ ' jects are between 50 and 80% of the total. A few galaxies, although characterized by a 



3 



compact morphology, are best fitted by an exponential distribution, more typical of local 
spirals. Our data also suggest that irregular EROs are found predominantly in the field, 
and that - on average - they tend to be characterized by the reddest colors. Finally, 
we plot the rest-frame Kormendy Relation (/i e vs. R e ) for a sample of 6 EROs with 
spectroscopic redshifts (z ~ 1.3), and estimate its evolution with respect to the local 
relation. 



1. Introduction 



Among the variety of objects discovered so far at high-redshift, a special class is rep- 
resented by the so called Extremely Red Objects (EROs hereafter), characterized by 
moderately faint near-IR magnitudes (K ^ 18 — 20), and extremely red optical infrared 



colors (e.g, R — K > 5, see for example Elston et al. |1988| , McCarthy et al. |1992| , Hu & 



Ridgway 1994 ). The observed colors and luminosities place this class of objects at z > 1, 
an hypothesis confirmed in a few cases by a direct spectroscopic measurement of the red- 



shift (Graham & Dey 1996, Spinrad et al. 1997, Stanford et al 1997 - S97 hereafter, Liu 
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can be either high-redshift starburst galaxies reddened by a large amount of dust, or 
passively evolving high-z ellipticals characterized by old stellar populations (> 1 Gyr). 

The importance of assessing the ERO nature and determining their space density is 
clear: the epoch of formation of massive elliptical galaxies is a crucial test for the standard 



hierarchical models for structure formation (e.g.: White & Recs 197S, Kauffmann et al. 



1993), which predict such objects to have formed relatively late from the merging of 



smaller-size objects (presumably disk galaxies). A large density of high redshift evolved 
ellipticals would imply severe revision to the hierarchical theories. The other relevant 
question is the global star formation history: calculations based on the observed rest 



frame UV flux (e.g. Madau et al. 1996 Connolly et al. 1997 ) might be significantly 
underestimated, if a large fraction of the overall star formation at high redshift takes 



place in highly obscured starburst galaxies (e.g. Steidel et al. 1999 ; Barger et al. 2000 ). 

One way to disentangle this ambiguity is provided, in some cases, by near-infrared 
spectroscopy, in particular if the ERO spectrum exhibits features revealing star- 
formation activity, such as the redshifted H a line; this is the case, for example, for the 
galaxy HR10 (Graham & Dey |1996| , Dey et al. |1999| ). More recently, deep near-infrared 
spectroscopy allowed to classify two more galaxies as likely starburst - although their 
spectra lack of spectral features - from the amount of reddening required to explain 
their overall spectral energy distribution (Cimatti et al. 1999] ). A different kind of test 
is provided by observations in the submm waveband, which traces the thermal emission 
by dust in the starbursts; this method was successful in the case of HR10 (Cimatti et 



al. 1998 , Dey et al. 1999 ) whose detection allowed its non-ambiguous classification as a 
dust reddened starburst, a result furtherly confirmed by the observation of its CO emis- 



sion (Andreani et al. 200C ). Other objects, first detected in the submm, have afterwards 



turned out to be EROs (Smail et al. |1999| Gear et al. |2000[ ). 

When images of sufficient spatial resolution are available, however, the most direct 
way to distinguish between the two classes is their morphology: elliptical galaxies are 
compact, regularly-shaped objects, whereas we expect starburst galaxy to look much 
more irregular (in particular, if the starburst is triggered by a merger, or if a large 
amount of dust irregularly distributed is present in the galaxy). HR10, imaged by HST, 



is consistently characterized by a clearly disturbed morphology (see Dey et al. 1999). 
For what concerns the total number density of EROs and their link with passively 



evolving ellipticals, the works by Cowie et al. |1994| and Hu & Ridgway |1994| suggest 
that at most a fraction of the present day ellipticals (~ 10%) could have its progenitors 



among EROs, but at present the question is far from being settled. Thompson et al. 1999 
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z ~ 2; Eisenhardt et al. 200C, finally, argue that the fraction of red galaxies at z > 1 might 
be significantly higher than previously thought, and consistent with a pure luminosity 



evolution scenario. Finally, Daddi et al. |200C| recently showed that EROs are strongly 
clustered and that such a clustering can explain the origin of the previous discrepant 
results on the surface density of z > 1 elliptical candidates as due to strong field-to-field 
variations. It has also been noted that even a small amount of star formation would 
drive a high-redshift elliptical galaxies towards bluer colors, so that it would be missed 
by a sample selection based on photometric properties only (for example, see Schade 



et al. 199S ). The problem of identifying high-redshift evolved galaxies, therefore, is not 
restricted to the ERO population alone; in this perspective, color-based selection criteria 



appear insufficient. Again, a different diagnostic tool (Franceschini et al. 1998, Schade et 



al. 1999) is provided by a quantitative analysis of morphological characteristics. A local 
elliptical galaxy is an evolved system from the point of view of both its stellar population 
and its internal dynamics; we may presume that, for some objects at high-z, a residual 
small star formation activity (and in general the overall stellar content of the galaxy) 
could affect the global colors but leave the shape of the brightness distribution more or 
less unchanged, so that such galaxies could be easily identified on a morphological basis. 
Of course this kind of approach requires imaging at high angular resolution, such as can 
only be obtained by space observatories (namely, by the Hubble Space Telescope - HST 
hereafter) . 



As a first effort to investigate the morphology of EROs, we present here a quan- 
titative analysis carried out on deep HST archive images both in the optical red and 
in the near infrared. Our aim is to identify elliptical galaxies using their morphological 
characteristics, and establish their fractional abundance with respect to the overall ERO 
population, assuming that their surface brightness distributions at z > 1 are similar to 
the ones observed in the local universe. In particular, at the resolution provided by HST, 
a first classification can be performed visually between compact and irregular objects; 
among the former ones, different distributions can then be distinguished by fitting differ- 
ent models to the data (for example, exponential and de Vaucouleurs profiles, typically 
associated to disk galaxies and ellipticals respectively). To this purpose, we have imple- 
mented a code for the analysis of the surface brightness distributions of such objects 
as observed by HST, tested its accuracy on a large number of simulated galaxies, and 
applied it to a sample of 41 EROs. 

This paper is organized as follows: we start describing our sample and the data avail- 
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resulting parameters and a morphological classification of the sample are carried out in 
Sect. [6] and ^; the conclusions follow in Sect. ||. 



2. The sample 

We selected a sample of candidate high-redshift ellipticals from a collection of EROs 
with published optical-infrared colors. We restricted our choice to the galaxies for which 
deep images in the red or near-infrared photometric bands were available from the HST 
archive - in particular observed by WFPC2 or NICMOS; roughly speaking, for objects 
at z > 1 these two instruments map respectively the UV and optical rest-frame spectral 
regions of the emitted radiation. Four of the sample galaxies are in the Hubble Deep 



Field South (HDFS), and were studied by Bemtez et al. [1999 



All the selected objects have I — K > 4 and/or R — K > 5; such thresholds are 
appropriate to select candidate elliptical galaxies at about z > 1, as shown in Fig. |], 
where we plot the theoretical evolution of the observer-frame colors, as a function of z, 
for typical stellar populations that we may expect to find in elliptical galaxies. 
The resulting set comprises 63 galaxies, but 22 of them are either hardly visible or not 
detected at all in the final reduced images (see the next section), and have therefore 
been excluded from the following analysis; these galaxies are listed in Table |[ The final 
sample of 41 EROs is listed instead in Table |[ together with some relevant information 
on each object; colors and K magnitudes from the literature are reported in Table ||. The 
-ftT-band magnitudes span rather homogeneously the range between 18 and 21, whereas 
the typical colors are 5 < R — K < 7 and 4 < I — K < 5.5. 

Clearly, the sample was not selected according to any fixed limit in total flux or 
surface brightness, but it rather comprises objects observed in different passbands and 
with different sensitivities^ as a consequence, it cannot be considerd complete at any 
flux level. On the other hand, the selection is based only on the availability of deep HST 
images, so that we do not expect any particular bias to be present; we also note that, 
in spite of its incompleteness, the size of this sample is unprecedented for this class of 
objects. Finally, since we are mainly interested in the structural characteristic of these 
objects, rather than in their intrinsic photometric properties, the lack of information 
about the redshift of most of the selected galaxies (9 spectroscopic and 6 photometric 
redshifts are available from the literature) does not represent a major problem. The 
surface brightness distribution of local elliptical galaxies exhibits little shape variation 



1 In Table |l| and |^ we report the limiting senstivity for the single frames, measured as the 
surface brightness (mag arcsec -2 ) at a 3-cr level on the single pixel; we estimated it as /un m = 
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Fig. 1. The curves show the R — K and I — K expected colors for passively evolving 
ellipticals. The models assume an exponential decaying star formation rate (SFR) with 
e-folding time r, solar metallicity, Salpeter IMF, and they are based on the 1997 release 



of the Bruzual & Chariot (1993) spectral synthesis models. The assumed cosmology is 



fi = 1 and -f/o=50 km s _1 Mpc -1 . The solid lines have r = 0.1 Gyr and Zf ormat i on =3,4,5. 
The dashed line has r = 0.3 Gyr and Zf ormat ion =5. The dotted lines at R — K = 5 and 
/ — K = 4 show the color thresholds that can be defined to select elliptical candidates at 
z > 1. Redder color thresholds allow the selection of higher-z ellipticals (e.g. R — K > 6 
or / — K > 4.5 correspond to z > 1.3). 



from the near ultraviolet to the near infrared so that similar passively evolving galaxies 
at z = 1 ~ 1.5 can be easily identified with our kind of analysis and with the photometric 
bands available. The possible effects of the different wavelength coverage of the WFPC2 
and NTCMOS images are discussed in section fi.3. 
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Table 1. The final sample 



ID 


R.A. (J2000) Dec. 


Ref. 


Field 


Name 


Instrum. 


Filter 


Exp. t. 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


01 


03:57:29.99 


+01:09:57.2 


[1] 


SSA4 


(28,19) 


WFPC2 


F814W 


12700 


24.8 


02 


06:50:51.98 


+41:30:39.3 


[2] 


4C41.17 


#8 


WFPC2 


F702W 


21600 


24.5 


03 


06:50:53.65 


+41:30:15.6 


[2] 


4C41.17 


#16 


WFPC2 


F702W 


21600 


25.4 


0-1 


08:48:30.79 


+44:53:34.8 


[3] 


GIG J 0848+44oo 


#237 


InIUVIUo 


T7 1 1 £ HIT 7 

r 160 VV 


11200 


25.0 


05 


08:48:32.42 


+44:53:35.1 


[3] 


CIGJ0848+4453 


#65 


NICMOS 


F160W 


11200 


25.0 


06 


08:48:32.99 


+44:53:46.6 


[3] 


CIGJ0848+4453 


#142 


NICMOS 


F160W 


11200 


25.0 


07 


08:48:35.99 


+44:53:36.1 


[3] 


GIG J 0848+44o3 


#70 




r 160 VV 


11200 


25.0 


08 


08:48:36.16 


+44:54:17.3 


[3] 


CIGJ0848+4453 


#108 


JNIGMOS 


F160W 


11200 


25.0 


09 


08:48:36.24 


+44:53:55.4 


[3] 


CIGJ0848+4453 


#135 


NICMOS 


F160W 


11200 


25.0 


10 


09:05:30.55 


+34:08:09.1 


[4] 


B2 0902+34 


Galaxy A 


WFPC2 


F814W 


5000 


23.8 


11 


10:19:32.85 


+05:34:35.8 


[5] 


MG 1019+0535 


Object D 


VV r PL* 2 


r 814 VV 


4700 


24.0 


12 


12:15:45.54 


—00:35:11.7 


[6] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


13 


12:15:46.57 


— 00:35:02.2 


[6] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


14 


12:15:47.38 


—00:34:45.2 


[6] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


15 


12:15:50.82 


-00:34:31.4 


[6] 


QSO1213-0017 


R7 


WFPC2 


F814W 


4500 


24.1 


16 


12:15:51.41 


—00:34:31.4 


[6] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


17 


12:15:53.10 


-00:34:43.8 


[6] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


18 


_LZ. 10. OZ. 1 1 


— UU. 00.114.0 


[6] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


19 


12:15:52.16 


-00:33:58.7 


T61 




R10 


V \ L JT \~>Z 




4500 


24 1 


20 


13:12:03.99 


+42:44:41.5 


[7] 


SSA13 


#144 


WFPC2 


F814W 


15300 


24.9 


21 


13:12:03.90 


+42:44:45.1 


[7] 


SSA13 


#143 


WFPC2 


F814W 


15300 


24.9 


22 


13:12:13.69 


+42:44:16.6 


[7] 


SSA13 


#79 


WFPC2 


F814W 


15600 


25.0 


23 


13:12:13.10 


+42:44:41.5 


[1] 


SSA13 


(34,2) 


WFPC2 


F814W 


15600 


23.8 


24 


13:12:14.11 


+42:43:55.9 


[7] 


SSA13 


#77 


WFPC2 


F814W 


15600 


25.0 


25 


13:12:14.49 


+42:45:15.0 


[1] 


SSA13 


K10 


WFPC2 


F814W 


15600 


25.0 


26 


13:12:16.74 


+42:45:08.7 


[1] 


SSA13 


K8 


WFPC2 


F814W 


15600 


25.0 


27 


13:12:15.91 


+42:44:43.0 


[1] 


SSA13 


(4,3) 


WFPC2 


F814W 


15600 


25.0 


28 


13:12:24.15 


+43:43:57.9 


[7] 


SSA13 


#138 


WFPC2 


F814W 


15300 


25.0 


29 


16:44:57.06 


+46:26:01.8 


[8] 


PC1643+4631 


HR14 


WFPC2 


F814W 


5300 


23.9 


30 


16:45:02.38 


+46:26:25.1 


[8] 


PC1643+4631 


HR10 


WFPC2 


F814W 


5300 


23.9 


31 


20:27:59.43 


-21:40:48.9 


[9] 


MRC2025-218 


Object a 


NICMOS 


F160W 


10259 


23.0 


32 


20:28:00.12 


-21:40:57.4 


[9] 


MRC2025-218 


Object c 


NICMOS 


F160W 


10259 


23.0 


33 


22:17:32.90 


+00:14:47.0 


[7] 


SSA22 


#107 


WFPC2 


F814W 


24000 


24.7 


34 


22:17:33.25 


+00:14:26.9 


[1] 


SSA22 


(23,-35) 


WFPC2 


F814W 


24000 


24.7 


35 


22:17:33.41 


+00:15:04.1 


[1] 


SSA22 


(21,0) 


WFPC2 


F814W 


24000 


24.7 


36 


22:17:33.75 


+00:15:29.8 


[1] 


SSA22 


(-33,27) 


WFPC2 


F814W 


24000 


24.7 


37 


22:17:41.06 


+00:15:01.2 


[1] 


SSA22 


#34 


WFPC2 


F814W 


14800 


24.8 


38 


22:32:50.70 


-60:38:30.3 


[10] 


HDFS 


ET4 


NICMOS 


F160W 


128441 


25.6 


39 


22:32:51.10 


-60:39:09.8 


[10] 


IIDFS 


ET3 


NICMOS 


F160W 


128441 


25.6 


40 


22:32:52.82 


-60:39:11.2 


[10] 


HDFS 


ET2 


NICMOS 


F160W 


128441 


25.6 


41 


22:32:55.46 


-60:38:32.4 


[10] 


IIDFS 


ET1 


NICMOS 


F160W 


128441 


25.6 



Notes to Table 1. 

Col. 4: references for the single objects. [1] Cowie et al. 1994 ; [2] Graham et al. 1994; [3] Stanford et al. 1997; [4] Eisenliardt & 
Dickinson 1992; [5] Dey et al. 1995; [6] Liu et al. 2000; [7] Cowie et al. 1996; [8] Hu & Ridgeway 1994; [9] McCarthy et al. 1992, 
Cimatti et al. 1999; [10] Bem'tez et al. 1999. Col. 6: Object name in the reference of Col. 4, or the coordinates in the published 
finding chart. Col. 7: The plate scale for the WFPC2 data is 0.0996 arcsec/pix, except for objects 2 and 23 (0.0455); for objects 
from 4 to 9 the scale is 0.2, for objects 31 and 32 it is 0.075, as well as for objects from 38 to 40. Col. 9: Exposure time in 
seconds. Col. 10: limiting surface brightness (AB mag arcesc -2 ), computed for the single pixel at a 3-<r level. 



3. Data reduction and analysis 

We have started our analysis from the pipeline-reduced HST datasets retrieved directly 
from the archive. Typically, for each field, the observations consist of a sequence of frames 
slightly shifted with respect to each other (shifts range from a few pixels to a few tenths), 
often strongly affected (in particular the WFPC2 ones) by cosmic ray hits. To obtain a 
single frame per object and to remove at the same time cosmic rays and residual bad 
pixels, we have performed a few further reduction steps to align and combine the images 
available for each field; a few tests carried out on field stars have shown that the procedure 
adopted to this purpose, outlined below, does not significantly degrade the quality of the 
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ID R.A. (J2000) Dec. Ref. Field Name Instrum. Filter Exp. t. //, im 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



42 


03:57:34.40 


+01:10:06.5 


[1] 


SSA4 


(-35,28) 


WFPC2 


F814W 


12700 


23.8 


43 


12:15:46.86 


-00:34:45.4 


[2] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


44 


12:15:47.43 


-00:34:38.4 


[2] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


45 


12:15:48.04 


-00:34:35.6 


[2] 


QSO1213-0017 


R6 


WFPC2 


F814W 


4500 


21.1 


46 


12:15:48.13 


-00:34:51.9 


[2] 


QSO1213-0017 




WFPC2 


F814W 


4500 


21.1 


47 


12:15:50.88 


-00:34:00.7 


[2] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


18 


12:15:51.56 


-00:34:13.8 


[2] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


19 


12:15:51.86 


-00:34:20.6 


[2] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


50 


12:15:51.87 


-00:34:54.1 


[2] 


QSO1213-0017 




WFPC2 


F814W 


4500 


24.1 


51 


21:07:13.07 


+23:31:19.6 


[3] 


4C23.56 


#14 


WFPC2 


F675W+Pol. 


6600 




52 


21:07:13.18 


+23:31:12.5 


[3] 


4C23.56 


#16 


WFPC2 


F675W+Pol. 


6600 




53 


21:07:13.95 


+23:32:08.0 


[3] 


4C23.56 


#35 


WFPC2 


F675W+Pol. 


6600 




54 


21:07:15.29 


+23:31:20.2 


[3] 


4C23.56 


#68 


WFPC2 


F675W+Pol. 


6600 




55 


21:07:15.46 


+23:31:37.8 


[3] 


4C23.56 


#70 


WFPC2 


F675W+Pol. 


6600 




56 


21:07:15.99 


+23:31:12.7 


[3] 


4C23.56 


#86 


WFPC2 


F675W+Pol. 


6600 




57 


21:07:16.09 


+23:31:23.1 


[3] 


4C23.56 


#92 


WFPC2 


F675W+Pol. 


6600 




58 


21:07:16.27 


+23:32:04.3 


[3] 


4C23.56 


#98 


WFPC2 


F675W+Pol. 


6600 




59 


21:07:16.44 


+23:31:12.6 


[3] 


4C23.56 


#102 


WFPC2 


F675W+Pol. 


6600 




60 


21:07:16.48 


+23:31:45.9 


[3] 


4C23.56 


#105 


WFPC2 


F675W+Pol. 


6600 




61 


13:12:19.36 


+42:45:09.0 


[1] 


SSA13 


(-33,27) 


WFPC2 


F814W 


15600 


25.0 


62 


13:12:17.58 


+42:45:04.9 


[1] 


SSA13 


(-14,23) 


WFPC2 


F814W 


15600 


25.0 


63 


22:17:33.06 


+00:14:59.8 


[1] 


SSA22 


(26,-4) 


WFPC2 


F814W 


24000 


24.7 



Notes to Table 2. All fields are at the scale of the Wide Field camera, except the one pertaining to object n. 42, imaged by the 
Planetary Camera. Columns are as in Table 1. Col. 4: references for the single objects. [1] Cowie et al. 1994 ; [2] Liu et al. 2000; 
[3] Knopp k Chambers 1997. Col 10: we do not compute a limiting surface brightness for the galaxies from reference [3], since 
they were observed in a filter+polarimeter configuration. 



The relative shifts between two different exposures of the same field have been eval- 
uated measuring the position of the peak in the cross-correlation of the two images. We 
have subsequently aligned all the frames using the DRIZZLE task (Fruchter & Hook 



199? ), preserving the original scale in the output images. We have chosen not to resam- 



ple the data on a smaller scale, since this yields an actual improvement of the spatial 
resolution in the final combined image only if the frames to be combined are shifted by 
non integer amounts, and if such shifts sample homogeneously the sub-pixel scale, which 
usually does not happen for our data. Also, as we will explain in the next section, we 
compare our data with model distributions convolved with a theoretical PSF, and the 
effects of the resampling on the PSF are rather difficult to quantify. In the end, a better 
accuracy in the evaluation of the PSF shape more than compensates a possible little loss 
in spatial resolution. The shifted frames are combined in a following stage, rather than 
by DRIZZLE itself, to achieve a more efficient cosmic-ray rejection. In the case of the 
HDFS, the public F160W image has been used, without any further processing. 

Using the ELLIPSE task in IRAF, a radial surface brightness profile has been ex- 
tracted for all the galaxies except one very irregular object (number 21 in Table 0). The 
final images of the sample galaxies are shown in Fig. || on a logarithmic scale, together 
with the respective radial brightness profiles. The map of object 21 is shown in Fig. pf we 
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Fig. 2. For each galaxy we show the final image in a logarithmic scale, and a plot with 
the radial surface brightness profile. The orientation of each map is with north up and 
east to the left, and the scale on the two axes is in arcsec; the filter is specified in the 
plot panel. Here, the dots and the solid line represent respectively the surface brightness 
profile of the galaxy and of the best fit model; both profiles are computed as the average 
intensity of the distributions along the same set of elliptical contours. For objects n. 5 
and 30 a satisfactory fit could not be obtained. Magnitudes are computed in the AB 
system. 
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4. Fitting the surface brightness distributions 

A model brightness distribution has been fitted to each object, using a modified version 



of the code described in Moriondo et al. 1998, Such a code was originally implemented 



to analyze the brightness distribution of nearby spiral galaxies by fitting to the data a 
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Fig. 2. Continued. 



4-1. The Point Spread Function 

The main changes introduced to make the fitting code suitable for the WFPC2 and NIC- 
MOS data concern the convolution of the model galaxy with the PSF, which is not axisim- 
metric and subject to significant changes from one point to another in the field. In other 
words, the gaussian approximation is no longer accurate enough, and the PSF needs to 
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Fig. 2. Continued. 

Such a theoretical PSFs proved to be accurate within about 15%, when compared with 
observed stars and considering the average residuals in the inner 5x5 pixels. The ac- 
curacy achieved by TinyTim is not worse than what could be obtained using stars in 
the field; this is mainly due to the variation of the PSF shape across the field of view 
(differences up to 20%, even in the central pixels, are easily observed), coupled to the 
fact that it is usually difficult for WFPC2 images - almost impossible for the NICMOS 
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Fig. 3. We identify object n. 21 with the one in the center of the map (the filter is 
F814W). We did not obtain a fit for this galaxy due to its disturbed morphology, and to 
the closeness of a bright companion. 

accuracy of the model convolution. This is particularly true in the case of WFPC2 data, 
where the image scale undersamples the PSF: a few tests have shown that in this case 
a good convolution of the model can be obtained only if it is performed on a finer grid. 
The convolved model is then rebinned to the proper scale and compared to the data. In 
the case of the HDFS, since the reduced frames have been resampled by DRIZZLE, we 
have chosen to use a PSF derived directly from the stars in the field. 

The convolution is performed as the inverse Fourier transform of the product of the 
direct transforms of the model and the PSF, using a Fast Fourier Transform algorithm. 



4.2. The model 

Because we want to identify elliptical galaxies, and given the small size of the sample 
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distribution in every fit is a generalized exponential (Sersic 1982| ): fi ~ exp(a„i? 1 / n ), 
including the case of an exponential distribution (n = 1) and of a de Vaucouleurs one 
(n = 4). The parameters of each fit are the effective radius R e and the effective surface 
brightness \i e of the distribution, as well as its center coordinates; the apparent ellipticity 
and position angle are held fixed since they can be determined more reliably from the 
ellipse fitting routine. The "shape index" n is also fixed, in every fit, to an integer value 
ranging from 1 to 6, due to the fact that, at the low signal-to-noise ratio (S/N) typical 
of our data, the fitting routine is not able to obtain a reliable estimate for both n and the 
other parameters at the same time. The best value of the shape index n for each galaxy 
is determined instead a posteriori, by choosing the least \ 2 resulting from the different 
fits. The accuracy of the final best-fit parameters, including n, has been assesed using a 
large set of simulated galaxies, and will be discussed in the next section. 

The fits are performed inside a circular region centered on the galaxy; its radius is 
chosen, using the radial brightness profile, as the one at S/N = 1 for the ellipse-averaged 
intensity. The background level is estimated on blank sky regions close to the source; its 
uncertainty turns out to be dominated, in most cases, by fluctuations on scales of order 
10 pixels or more, due to a non-perfect image flattening. 



5. Simulations of faint galaxies 

To establish the accuracy that can be attained by the fits for all the relevant parameters, 
we have tested our code on a large set of simulated galaxies. We have chosen model 
distributions spanning the same range in effective radius and total flux - in terms of 
instrumental units - as our sample galaxies, and a range of values (from 1 to 5) for the 
shape index. Since most of the objects considered have been imaged with the wide field 
camera of WFPC2, we have convolved the test distributions with a typical wide-field 
PSF. In this case, the plate-scale is the one which mostly undersamples the PSF itself, 
yielding the worst conditions for the retrieval of the correct parameters' value. We expect 
therefore the uncertainties derived from our simulations to be basically correct for the 
WFPC2, and conservative estimates for the other types of data: a few simulations carried 
out with the characteristic sampling of the Planetary Camera and of NICMOS show that 
this is indeed the case. 

Noise has been added to the simulations at a level typical for our data, both on the 
pixel scale and at smaller spatial frequencies. The resulting images have been finally fitted 
using our code, allowing for some error in the estimate of the PSF and using different 
trial values for n, as in the case of the real galaxies. 
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each value of n, from 1 to 4. The dotted contours map the S/N, evaluated theoretically, 
considering the total flux inside the isophotal radius at S/N = 1 on the single pixel. 
These values do not account for the effect of the PSF, and are therefore less reliable at 
low surface brightness levels and for sizes of the order of the PSF width (say, log(i? e ) ~ 0). 
For this reason, the detection limits for our data have been deduced a posteriori from 
the simulations, and are represented by the dashed lines in the lower right corners of the 
plots. Such limits appear to depend on n; this is partly explained by the fact that, at fixed 
effective radius and surface brightness, there is a slight increase of S/N with increasing n. 
However, the greatest contribution (about 80%) is due to the steeper central peak that 
characterizes the large-n distributions, making them more visible in the background 
noise. 



5.1. Disentangling the n — 1 and n — 4 distributions 

A first check can be carried out assuming that all elliptical galaxies are characterized 
by a de Vaucouleurs brightness distribution (n = 4); we can consider this as a first 
order classification, since nearby ellipticals show in fact a variety of shapes, that can be 
quantified by different values of the exponential index n ranging from about 2 to 10 and 
higher (see for example Caon et al. 1993| , or Khosroshahi et al. |2000 ). If we restrict our 
test to the simulations with n — 1 and n — 4 only, we find that the correct value of n 
can be retrieved, on the basis of a y 2 estimator, for the whole parameter space explored. 
Such a simplified classification, therefore, is possible for all of our objects. 



5.2. Systematic errors 

In the following step we assigned to every galaxy its best-fit n value, choosing from 1 to 
6, and compare the derived parameters with the original ones. We can start our analysis 
of the results by looking for systematic trends. Starting with our estimates for the index 
n, if we plot the average measured values vs. the true ones (Fig. |E|), we actually observe 
a tendence for the higher n's to be underestimated; in particular we find 



n true = (1.37 ±0.30)n trae + (-0.41 ±0.35) . (1) 

This effect depends only slightly on the choice of the PSF; it is probably due to the fact 
that, increasing n at fixed total flux, the peak of the distribution gets sharper, and its 
wings fainter and wider, but these differences tend to be concealed by the effect of the 
PSF on one hand, and by the backgroud noise on the other. 
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Fig. 4. The region covered by the simulations, and the estimated S/N ratio. The different 
panels show the results for different values of n, from 1 to 4. Each panel represents the 
R e -fi e plane in instrumental units (pixels and counts per pixel). The region covered by 
our simulations, including the area occupied by the sample EROs, lies above the solid 
line; its slope is at constant total flux. The level marked by the dotted horizontal segment 
is the noise level of the background at 1 a, whereas the dashed lines in the lower right 
corner represent the detection limits for the n — 1, 2, 3, 4 distributions in the respective 
panels. The dotted contours map the estimated S/N at the values indicated by the labels. 



level) tend to be slightly overestimated, and large values tend to be underestimated. The 
behaviour is very similar for all the n values, so that we can adopt average corrections: 



log(iW) - (1-15 ± 0.05) log(i? obs ) + (-0.09 ± 0.03) 



(2) 
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2 4 
Best n (real) 

Fig. 5. The trend of measured vs. real n, for the simulated galaxies, n values larger than 
2 tend to be underestimated. 

Since these latter corrections are significant at a 3 er level, whereas Eq. [I] is significant 
only at 1 a, and since applying Eq. [I] to the estimated n's would lead to non-integer 
values for this parameter, we choose to correct only R e and /i e and leave n unchanged, 
keeping in mind that n values greater than two might be somewhat underestimated. 

5.3. Mapping the parameter space 

We turn now to examine how accurately the relevant parameters are retrieved in the 
various regions of the parameters' space, starting with the shape index n. Fig. ^ shows 
again the R e -fJ. e plane; in this plot the dots in each panel represent the estimated location 
of our simulated galaxies for the different values of n. The accuracy with which n can be 




Fig. 6. Accuracy in the estimates of n. The plot is similar to the one in Fig. |j. In 
each panel the dots are placed at the estimated location of the simulated galaxies, with 
their size quantifying the accuracy in the estimate of n: the small dots correspond to the 
correct value, the medium-size ones imply an error of ±1, the large ones an error grater 
then 1. In the region above the dotted line, we can relibly distinguish between n = 1 and 
n > 1 distributions. 

We find that the correct value of n is retrieved in most cases for the n = 1 and n = 2 
models; the error for the n = 3 and n = 4 ones is more typically 1 in large portions of 
the plane, partly due to the systematic effect described previously. As a consequence, 
exponentials are almost always recognized as such, so that if the best fit is for n ^ 1, the 
distribution is certainly non-exponential. As expected, for all values of n, low flux and 
low surface brightness objects tend to be affected by larger errors. The main conclusion, 



18 Moriondo et al.: ERO morphology 

is the locus where both exponentials and n = 2 distributions are recognized as such, and 
larger n distributions are affected at most by an error of 1. A comparison with Fig. |^ shows 
that the limit roughly spans S/N values between 10 and 80. We have checked this result 



using the theoretical approach described in Avni 1976 when one or more parameters 
are evaluated via a x 2 minimization, the method allows to assign a confidence level to 
each parameter relying on the variations of the \ 2 around the minimum in the parameter 
space. Although the computations are exact only in the case of linear fits, the method 
provides anyway a useful check on our findings; indeed, we find that our estimates for 
the uncertainty of n are broadly consistent with the ones evaluated theoretically for a 
90% confidence level. In particular, the Avni method confirms that in the area above the 
dotted line, exponential and non-exponential distributions can be reliably distinguished. 

A mapping of the parameter space, analogous to the one plotted in Fig. [| has been 
produced also to estimate the uncertainties on R e and fx e , corrected according to Eqs. 
|^ and ||. An intersting result is that the derived errors are relatively independent of the 
estimate of n, in the sense that a wrong estimate of the shape index does not necessarily 
mean larger errors for R e and fi e . Most likely, whereas the choice of n is influenced mainly 
by the accuracy of the PSF, the estimates of R e and ^ e are more strictly related to the 
quality of the fit a whole; in other words, if a wrong n may compensate for the effect 
of a wrong PSF, a good estimate for R e and /i e can be achieved anyway, as long as the 
quality of the fit is good. 

For what concerns the values of the ellipticity and position angle (that are fixed a 
priori), we find that the typical errors associated to their estimates do not affect signif- 
icantly the accuracy of the output parameters (center coordinates, effective radius and 
surface brightness), nor the choice of the best n value. We estimate the typical errors 
on the ellipticity to be around 0.1, and from 5 to 10 degrees for the position angles. 
The center coordinates are usually determined with great accuracy (~ 0.3 pixels): due 
to the asymmetries in the psf, this is better than what can be achieved by fitting an 
ellipse-averaged profile to the galaxies. 

To summarize, we have tested our fitting code on a large set of simulated galaxies, 
and assessed the accuracy that can be attained for the various relevant parameters in the 
region of the parameters' space covered by the real data. In the next section, the results 
presented so far will be used to estimate the errors on the parameters derived for each 
galaxy. 
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6.1. Irregular and compact objects 

The first classification that can be carried out for our sample is between compact, isolated 
objects, and galaxies that are clearly undergoing an interaction or exhibit an irregular, 
diffuse shape. The easiest way to do this is by visual inspection, since all our objects 
are well resolved. We find 35 compact galaxies out of 41, or 85% of the sample. The 



ERO HR10 (Hu & Ridgway |1994| ; Graham & Dey |1996| ) is included in the subsample 
of irregular objects. We have attempted to recognize in each of these irregular objects a 
brighter component, such as could be expected in a merging system, and obtain a fit of 
its brightness distribution after a proper masking of the surrounding areas. This was not 
possible for object 21, object 5, and object 30 (HR 10): in the case of object 21, as we 
mentioned in Sect. |^, the surface brightness distribution is too diffuse to isolate a major 
component; object 5 is compact, but its nucleus has a rather irregular shape, probably 
due to the superposition of two or more close and equally bright components; in object 
30, a brighter component can be easily recognized, but we did not obtain a satisfactory 
fit to its brightness distribution. A radial brightness profile was anyway extracted for 
object 5 and object 30, as shown in Fig. ||. A few more galaxies are close to other objects 
that, however, are not disturbing their morphology: in these cases a proper masking was 
also applied, as indicated in Table ^, to avoid any problem with the fitting procedure. 



6. 2. Structural parameters 

The best-fit profiles are plotted in Fig. |[ the best-fit parameters for each object are shown 
in Table ||, while Fig. (?] shows the estimated location in the R e -^e plane of all our sample 
galaxies with different symbols for each instrumental configuration. The encircled symbols 
correspond to the 3 galaxies classified as irregular/interacting for which we could obtain 
a fit, after masking the lower flux companion: since for these objects we consider only a 
fraction of the total flux, they are typically placed in the lower part of the plot. Effective 
radii and surface brightnesses are plotted in instrumental units (pixels and counts/pixel), 
to allow a direct comparison with Fig. |^. The y coordinates of the data points, however, 
are not exactly the values determined by the fit, since each galaxy has been scaled to the 
noise level adopted for the simulations by applying a proper shift along the brightness 
axis. The uncertainties listed in Table [| are evaluated by interpolating the results from 
the simulations at the locations of the real galaxies in the parameters' space; we report 
1-tr errors for all the parameters except n: when its integer value is retrived correctly, 
we assign to this quantity a formal error of 0.5, otherwise the integer value reported 
corresponds to the largest possible error. As we mentioned in the previous section, we 
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Table 3. Photometric and structural parameters 



ID n R e pi e e M ZP R - K I - K K Notes 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



01 


2 


(2 


0) 





16 


(0.01) 


23.81 


(0.59) 





7 


26.13 


(0.06) 


22 


10 




1.9 


20.3 


C 


02 


5 


(1 


0) 


1 


07 


(0.29) 


26.55 


(0.53) 





2 


23.14 


(0.23) 


22 


72 


6.16 




18.41 


C 


03 


3 


(1 


0) 





29 


(0.03) 


24.44 


(0.03) 





2 


24.16 


(0.02) 


22 


70 


6.16 




18.41 


C 


04 


1 


(0 


5) 





29 


(0.02) 


21.37 


(0.10) 





3 


21.72 


(0.01) 


22 


85 


5.19 




19.53 


CM 


05 






























5.86 




18.11 


I 


06 


2 


(0 


5) 





14 


(0.02) 


19.59 


(0.02) 





3 


21.16 


(0.02) 


22 


85 


5.91 




18.90 


c 


07 


4 


(1 


0) 


1 


10 


(0.13) 


22.96 


(0.28) 





2 


19.61 


(0.10) 


22 


85 


6.09 




18.14 


c 


08 


5 


(1 


0) 





79 


(0.08) 


22.90 


(0.25) 





2 


20.15 


(0.08) 


22 


85 


5.87 




18.47 


CM 


09 


4 


(1 


0) 





68 


(0.07) 


23.09 


(0.20) 





2 


20.71 


(0.08) 


22 


85 


6.21 




18.86 


c 


10 


1 


(0 


5) 





35 


(0.02) 


23.54 


(0.21) 





1 


23.65 


(0.04) 


22 


10 


6.1 




18.3 


CM 


11 


4 


(3 


0) 





74 


(0.12) 


26.47 


(0.02) 





3 


24.11 


(0.02) 


22 


10 


> 7 






IM 


12 


2 


(2 


0) 





11 


(0.01) 


22.50 


(0.04) 





4 


24.76 


(0.03) 


22 


08 


> 6 






c 


13 


2 


(0 


5) 





20 


(0.02) 


23.20 


(0.25) 





5 


24.36 


(0.01) 


22 


08 


5-6 






c 


14 


3 


(1 


0) 





29 


(0.04) 


23.34 


(0.21) 





3 


23.18 


(0.13) 


22 


08 


5-6 






c 


15 


2 


(0 


5) 





36 


(0.02) 


24.33 


(0.32) 





1 


23.64 


(0.04) 


22 


10 


5.84 




18.04 


c 


16 


1 


(1 


0) 





10 


(0.02) 


22.56 


(0.32) 





1 


24.99 


(0.01) 


22 


10 


5—6 






c 


17 


3 


(1 


0) 





33 


(0.01) 


22.40 


(0.23) 





3 


21.93 


(0.02) 


22 


10 


5—6 






c 


18 


1 


(0 


5) 





24 


(0.02) 


23.37 


(0.21) 





3 


24.19 


(0.01) 


22 


10 


5-6 






c 


19 


4 


(1 


0) 





39 


(0.06) 


24.06 


(0.21) 





1 


23.25 


(0.14) 


22 


10 


6.26 




18.19 


c 


20 


2 


(0 


5) 





10 


(0.01) 


22.15 


(0.41) 





2 


24.39 


(0.06) 


22 


10 




4.7 


19.5 


c 


21 
































4.7 


19.3 


I 


22 


1 


(0 


5) 





45 


(0.03) 


24.12 


(0.25) 





6 


24.14 


(0.05) 


22 


01 




6.2 


19.3 


c 


23 


2 


(2 


0) 





10 


(0.02) 


23.88 


(0.47) 





3 


26.28 


(0.01) 


22 


0,5 




4.9 


20.6 


c 


24 


2 


(0 


5) 





11 


(0.01) 


22.51 


(0.01) 





4 


24.77 


(0.01) 


22 


01 




4.7 


20.0 


c 


25 


1 


(1 


0) 





14 


(0.02) 


23.85 


(0.39) 





2 


25.72 


(0.02) 


22 


04 




5.1 


19.9 


c 


26 


1 


(0 


5) 





25 


(0.04) 


25.51 


(0.16) 





1 


25.91 


(0.16) 


22 


01 




4.9 


19.7 


CM 


27 


1 


(1 


0) 





22 


(0.04) 


25.86 


(0.07) 





2 


26.65 


(0.05) 


22 


14 




4.8 


20.7 


c 


28 


1 


(1 


0) 





60 


(0.10) 


26.28 


(0.02) 





2 


24.24 


(0.02) 


22 


10 




4.8 


19.4 


c 


29 


1 


(2 


0) 





14 


(0.02) 


23.61 


(0.37) 





2 


25.46 


(0.02) 


22 


09 




6.2 


18.7 


c 


30 
































6.5 


18.42 


I 


31 


5 


(1 


0) 





51 


(0.09) 


21.53 


(0.09) 





3 


19.89 


(0.06) 


23 


11 


5.1 




17.9 


c 


32 


3 


(1 


0) 





53 


(0.08) 


21.68 


(0.02) 





2 


20.05 


(0.02) 


23 


11 


5.1 




18.1 


c 


33 


5 


(1 


0) 





42 


(0.11) 


24.82 


(0.27) 





2 


23.45 


(0.25) 


22 


14 




4.5 


18.8 


c 


34 


1 


(0 


5) 





40 


(0.02) 


25.00 


(0.27) 





3 


24.70 


(0.04) 


22 


11 




7.7 


20.3 


IM 


35 


3 


[2 


0) 





11 


(0.03) 


23.82 


(0.37) 





2 


25.46 


(0.22) 


22 


11 




4.5 


20.2 


c 


36 


1 


(1 


0) 





18 


(0.03) 


24.91 


(0.45) 





1 


26.07 


(0.12) 


22 


14 




5.1 


20.4 


c 


37 


2 


(2 


0) 





27 


(0.06) 


25.07 


(0.10) 





3 


25.24 


(0.06) 


22 


10 




4.6 


19.6 


IM 


38 


4 


(1 


0) 





47 


(0.10) 


24.60 


(0.33) 





3 


23.25 


(0.18) 


22 


85 


8.2 


7.6 


21.06 


c 


39 


3 


(1 


o) 





10 


(0.01) 


19.93 


(0.08) 





6 


21.81 


(0.03) 


22 


85 


6.7 


4.9 


19.55 


c 


40 


2 


(0 


5) 





28 


(0.02) 


24.23 


(0.33) 





5 


24.63 


(0.07) 


22 


85 


5.1 


4.3 


22.33 


CM 


U 


4 


1 


0) 





27 


(0.03) 


20.98 


(0.20) 





6 


21.37 


(0.08) 


22 


85 


5.3 


1.1 


18.97 


CM 



Notes to Table 3. 

Col. 2: best-fit shape index n. Col. 3: effective radius in arcsec. Col. 4: effective surface brightness in mag arcsec -2 (AB system). 
Col. 5: apparent ellipticity. Col. 6: total magnitude. Cols. 8, 9: published colors, in the Vega system. For objects 12-19, the 
reported colors are 1Z — K] following Steidel k. Hamilton 1993, and assuming R — I — 1.5, we have R = 1Z + 0.18. For objects 
31-32 we report the published colors from Cimatti et al. 1999 and Daddi et al. (in preparation), which are different from the ones 
previously measured by McCarthty et al. 1992. For objects 38-41, the reported colors are converted to the Vega photometric 
system from the estimates (in the AB system) by Lanzetta et al. 1998; we have assumed (R— K)ve 9 a = (R — F222M) ab + 1.48. 
Otherwise, the colors are extracted from the references in Table 1. Col 10: published K magnitudes, in the Vega system. No 
values are available for objects 11-14 and 16-18. For objects 29, 30 K' magnitudes are reported. For objects 38-41, the reported 
values are converted to the Vega photometric system from the estimates (in the AB system) by Lanzetta et al. 1998; we have 
assumed Kvega — Kab — 1-76. Col. 11: morphological classification. C — compact; I — irregular/interacting; M — masked. 



irregular) the thoretical estimate exceeds the one derived from the simulations; for these 
galaxies the uncertainties reported are the thoretical ones. 



The 4 HDFS galaxies (38-41) were studied by Bemtez et al. 1999 with a technique 
very similar to ours (a best fit to the brightness profiles with a de Vaucouleurs' law), so 
that their and our results can be easily compared. We find that a de Vaucouleurs law is 
the best fit to the data for two of these galaxies (n = 4 ± 1), the other two being best 
represented by an n — 2 ± 0.5 and n = 3 ± 1 profile respectively. For what concerns the 
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differences between the two works are not relevant, and characterized by only a modest 
scatter in the measured quantities. Excellent agreement is then found, in the case of 



object 39, with the results by Stiavelli et al. 1998: although our best fit for this galaxy is 



with n = 3, both our effective radius and our total flux are coincident with the Stiavelli 
et al. values, derived adopting a de Vaucouleurs' distribution. 

Figure ^ shows again the \x e -R e plane in standard units (arcsec, mag arcsec -2 ), with 
a different symbol for each filter; the dotted line represents the slope of constant flux, at 
fixed shape index n. The plot illustrates the limits in size and surface brightness of the 
sample in the HST filters. 

6. 3. The shape index n and the fraction of ellipticals 

We evaluated previously, through visual inspection, the fraction of irregular objects, 
concluding that they constitute only a minority of our ERO sample. For what concerns 
the shape of the best-fit distributions, we performed two types of classifications. 

A first order classification was performed by comparing the results assuming that 
each galaxy can be properly described by either an exponential distribution (n = 1) or 
a de Vaucouleurs one (n = 4). To do that, we considered only the simulations belonging 
to these two classes, as we have seen in the previous section that the true n value can 
be retrieved for all the galaxies in the sample. The resulting number of de Vaucouleurs 
distributions is then 21 out of 41 (51%). 

A more detailed classification was made leaving n free to vary among the integer 
values n — 1, 2, 3, 4, 5, 6. Using this approach, the relative abundance of non-exponentials 
{n > 2) is slightly higher if we choose the best n value for each galaxy from the whole 
set of fits performed: 25/41 (61%). In particular, four galaxies previously catalogued as 
exponentials are now fitted better by an n = 2 distribution, whereas the other 10 objects 
with n = 1 have their best fits confirmed. As we discussed previously, for the galaxies 
placed above the dotted line in Fig. ^ we can reliably distinguish between n — 1 and 
n =/= I, whereas at fainter fluxes n > 1 distributions may be mistaken for exponentials. 
This "high-signal" subsample, therefore, provides a particularly accurate estimate of the 
fraction of likely ellipticals which, in this case, is even larger, amounting to 81% (21 out 
of 26). 

As mentioned in Section 2, for z > 1, the WFPC2 and NICMOS images cover the 
rest-frame UV and the optical spectral regions respectively. As it is well known that the 
galaxy morphology depends strongly on the wavelength (e.g. Kuchinski et al. 2000] and 



references therein), one may argue if this can have effects on our results. In this respect, 




Fig. 7. The location of our sample galaxies in the fi e -R e plane; the units are counts 
per pixel and pixels respectively. Different symbols correspond to different instruments 
and/or detectors: filled circles for NICMOS camera 3; open circles for NICMOS camera 2; 
filled squares for the HDFS galaxies; open triangles for Wide Field data; filled triangles for 
Planetary camera data. The encircled points are those classified as irregular/interacting 
in Table ||. The other symbols are as in Fig. |[ 



as elliptical both in WFPC2 and in NICMOS images. Second, if a galaxy is irregular, 
then it would be reliably classified as such both in WFPC2 and in NICMOS images (e.g. 
HR10; Graham & Dey 1996[ Dey et al. 1999 ). Finally, there could be cases of elliptical 
galaxies with a disk component having n = 1 if observed in the optical (WFPC2) and 
n > 2 if observed in the near-IR (e.g. spheroidal galaxies with a disk component becoming 
more prominent in the rest-frame UV) . Our analysis does not allow to investigate if such 
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Fig. 8. The location of our sample galaxies in the fi e -R e plane; the units are mag 
arcsec -2 and arcsec respectively. In this plot different symbols for the data points corre- 
spond to different filters: filled circles for F160W, open triangles for F8I4W, open squares 
for F702W. The encircled points are those classified as irregular/interacting in Table ||. 
The slope of the dotted line is at constant flux. 



To summarize, we conclude that the galaxies that can be reliably classified as el- 
lipticals amount to 50-^80% of the total sample of 41 EROs. Although our sample is 
incomplete, we note that our results are in good agreement with those of Stiavelli & Trcu 



200C based on a complete sample of NICMOS-selected EROs. 

Such a high fraction of ellipticals strengthens the scenario proposed by Daddi et 
al. |000 who suggested that, because of their strong clustering, EROs are likely to be 
dominated by ellipticals rather than dusty starbursts. 
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Table 4. ERO morphology: a summary 



Whole sample 


Total 




Irr. 


Exp. (n = 1) 


Ell. (n = 4) 


All 


41 


6 


(15%) 


14 (34%) 


21 (51%) 


Cluster 


14 


1 


(7%) 


5 (36%) 


8 (57%) 


Field 


27 


5 


(19%) 


9 (33%) 


13 (48%) 


High signal 


Total 




Irr. 


Exp. (n = 1) 


Ell. (n > 2) 


All 


26 


1 


(4%) 


4 (15%) 


21 (81 %) 


Cluster 


11 





(0%) 


2 (18%) 


9 (82 %) 


Field 


15 


1 


(7%) 


2 (13%) 


12 (80 %) 



7.1. Field and cluster objects 

Since our sample includes galaxies both from the field and from a cluster environment, 
we can investigate the eventual differences between these two subsamples. Following the 
conclusions of the authors, we will consider as cluster members the objects previously 



studied by Liu et al. 2000 and S97, and assume the rest of the sample to be representative 
of the ERO field population; the two subsamples include 14 and 27 objects respectively. 

The most remarkable difference between them is that most of the galaxies classified as 
irregular (5 out of 6) belong to the field population. On the other hand, the only irregular 
cluster object (n. 5 in Table is neither diffuse nor characterized by two interacting 
components, but is rather a compact galaxy with an irregular core; also, its spectrum 
does not exhibit features typical of ongoing star formation (S97; again, one example of 
an apparently "old" object that does not resemble local ellipticals). We conclude that, if 
some starburst galaxies are present among the EROs, they are not likely to be found in 
clusters. Considering only the field population, the fraction of irregular galaxies is only 
slightly increased (19%) with respect to our previous estimate. 

For what concerns the fraction of non-exponential profiles, it is roughly the same for 
the two subsamples, again close to 80% for the high-signal objects. 

In Table |] we report a summary of the sample statistics. In the upper half of the 
table we consider the whole sample, divided into irregular galaxies, exponentials, and de 
Vaucouleurs. In the lower half we consider only the high-signal subsample, for which we 
distinguish between n = 1 and n > 1. 



7.2. Red exponential galaxies 
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ellipticals if we use the local objects as a reference; on the other hand, the regularity of 
their surface brightness distributions tends to exclude the hypothesis of heavily reddened 
objects. The existence of this subclass, therefore, implies that the ERO population is 
apparently more composite than previously thought, a conclusion also emerging from 



the work by Liu et al. |2000[ Stiavelli & Treu |2000| and Corbin et al. |2000| . The possibility 
that such objects are undergoing an intermediate post-merging phase that eventually 
ends up in an elliptical galaxy is in contrast with the simplest monolithic collapse model, 
in which all ellipticals are formed at high redshift. Its implications in the framework 
of the different scenarios for galaxy formation certainly deserve further, more detailed 
investigation. At the same time, we cannot exclude that some ongoing star formation, 
suitably distributed throughout the galaxy, might transform an elliptical-like bringhtness 
distribution into one of the kind observed. Again, this hypothesis could be tested by a 
better characterization of the stellar content of these exponential objects. 



7.3. Morphology and colors 

In Fig. ^ we plot the median I — K colors, computed for the different morphological 
classes introduced: irregular objects (corresponding to the n = point), exponential 
distributions, and the ones with n > 1. To include the galaxies with only R — K available, 
we adopted a color R — I = 1.5 to convert the R — K color into I — K . The adopted 
R — I = 1.5 is an average value estimated for elliptical galaxies using the synthetic 
spectral-energy distributions of Bruzual & Chariot (1997) for a set of passively evolving 
models at 1 < z < 2 (see the caption of Fig. 1). Since the R — I color depends strongly 
on where the 4000 A break falls, the adopted R — I = 1.5 has an associated uncertainty 
of ±0.3 mag. 

The dotted line represents the average I — K for n > 2 galaxies. Although each 
morphological bin is characterized by a significant scatter, there seems to be a tendence 
for the irregular objects to be characterized by the most extreme colors, with an average 
I — K exceeding the mean color of likely ellipticals (n > 2) by about 1 magnitude. 
Inverting the argument, we find that among the 6 reddest EROs (I — K > 5.5) only one 
exhibits the typical morphology of an elliptical galaxy. 

The adopted R — I = 1.5 is most likely appropriate for high redshift ellipticals, but 
might be somewhat overestimated for possible starburst galaxies; thus, if we adopt a 
smaller R — I for the irregular objects, the difference between them and the rest of the 
sample is furtherly increased. This result does not change also if we consider only the field 
subsample, and it is broadly consistent with the recent findings that the EROs detected 
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Fig. 9. Median I — K colors for the different morphological classes; the n = point 
corresponds to the irregular objects of the sample. The horizontal dotted line is the 
average I — K color for the n > 2 galaxies. Irregular galaxies are characterized, on 
average, by the reddest colors. 



detailed investigations - is therefore suggested between morphology, submm emission, 
and optical/infrared colors (in particular / — K): high z ellipticals and starburst might 
in fact exhibit a different behaviour with respect to each of the three parameters. 

7.4- The distribution of n 

We can consider the subsample of galaxies whose best shape index n is > 2 as a set of 
likely high-z ellipticals. For these objects we can compare the distribution observed for n 



with the results found for local samples of elliptical galaxies. Caon et al. 1993 determined 
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the range 2-8. More recently, the same kind of distribution was derived by Khosroshahi 



et al. 2000 for a sample of elliptical galaxies in the Coma cluster. 

In the top panel of Fig. ^ we plot our distribution for the high-signal subsample 
with no correction applied to the derived n values (thick solid line). The dotted line 
is the same distribution, corrected for the systematic effect described in Sect. || (Eq. 
|l|). In the bottom panel we plot the distributions for the Virgo cluster and the Coma 
cluster. Quite surprising, the two "local" histograms appear rather different, with the 
Virgo distribution extended up to large n values, and clearly peaked at n = 2, and the 
Coma distribution characterized by a peak and an upper cutoff at n = 4. Due to this 
diversity a comparison with our data is rather difficult; we just note that our "high- 
z" distribution looks somewhat intermediate between the two, being quite flat between 
n = 2 and n = 5, and confined to this range. Without attempting any deeper comparison, 
we limit ourselves to consider this result as consistent with our claim that we are actually 
looking at a population of elliptical galaxies. 

7.5. The Kormendy relation 

Scaling relations represent a powerful tool to investigate the evolution of galaxies at 
high rcdshift; in the case of elliptical galaxies, the Kormendy Relation between effective 
surface brightness and radius is relatively easy to build for a sample of ellipticals, when a 
detailed analysis of the brightness distributions and an accurate measure of the redshifts 



are available. Previous studies of this kind (for example Fasano et al. 199? , Ziegler et 



al. 1999 - Z99 hereafter, Roche et al. 1998]) observed, as expected, an increase of the 



rest-frame surface brightness with redshift, but the type of evolution implied (passive or 
partially active) is not yet well constrained by the models. 

For 6 of our n > 2 compact galaxies a spectroscopic measure of the redshift has 
been published (in particular, 4 galaxies from S97 and 2 from Liu et al. [2000 ). All these 



galaxies are likely to reside in a cluster environment, and all of them are approximately 
at z = 1.3. As a consequence, they make up a particularly homogeneus set, well suited 
to pinpoint a particular time in the luminosity evolution of cluster ellipticals. Four more 
elliptical candidates (the ones in the HDFS) have photometric redshifts measured, but 
the two estimates available (Chen et al. 1998| , Benitez et al. 1999) are quite discrepant 



and they have been excluded from the following analysis. 

We have used the 6 spectroscopic redshifts to derive the rest-frame parameters of 
the relative galaxies, following the prescriptions outlined in Z99; in particular, the ob- 
served F814W and H160W surface brightnesses have been corrected for the cosmological 
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Fig. 10. Top panel: The distribution of the shape index n derived from our sample, 
both uncorrected (thick solid line) and corrected (dotted line) for the systematic effect 
described in Sect. ||. Bottom panel: the distributions derived by Caon et al. 1993 (dotted 



line), and by Khosroshahi et al. 2000 (dashed line) 



models described in Pozzetti et al. 1998. The observed luminosities have been corrected 



for the Galactic extinction using the results by Schlegel et al. 1998 In Fig. |TT] we plot 
the rest-frame Kormendy Relation for the 6 selected galaxies, derived adopting Ho = 50, 
qo = 0.05, and a single stellar population model with redshift of formation Zf — 5. As a 
local reference, the solid line is the relation reported by Z99 for the whole sample studied 



by J0rgensen et al. [1994 



He = 3.46 (+1.17, -0.73)i? e + 19.46 (+0.33, -0.55) . (4) 
An upwards shift of the data points with respect to the local relation is evident: a 
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Fig. 11. The B-band, rest-frame Kormendy Relation for a subsample of "cluster" EROs 
at z <~ 1.3. The rest-frame parameters are evaluated assuming Ho — 50 km s _1 Mpc -1 
and go = 0.05. The solid line represents the local relation; the dotted line is the best fit 
to the data, with the slope fixed at the local value; the dashed line is a fit with both 
parameters free. The typical uncertainty on the data points is shown in the upper right 
corner of the plot. 

template (see Z99 for the details), or increasing the adopted value of q a up to 0.5. The 
best fit with both parameters free (dashed line) is 

fi e = (4.1 ± 0.5)i? e + (17.45 ± 0.15) , (5) 

consistent with the constant-slope hypothesis. The possible slight steepening of this 
relation should be considered with caution, both because of the very few data points, 
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A comparison with published results shows that the measured shift in surface bright- 
ness is consistent with the predictions of evolutionary models for elliptical galaxies (for 
example, with the Pure Luminosity Evolution models considered by Roche et al. 199§| ), 
as well as with the trends of luminosity and surface brightness vs. z observed at lower red- 
shifts. Z99, for example, find a difference of 0.8~0.9 i?-mag for two clusters at z = 0.55, 



whereas Schade et al. 1999 estimate a luminosity evolution of about 1 mag for a sample 
of field galaxies at redshift between 0.75 and 1. 



8. Summary and conclusions 

We have implemented a code to analyze the surface brightness distribution of faint galax- 
ies, optimized to work on deep HST images, and we have tested its performances on a 
large set of simulated galaxies. We have then used this code to fit model brightness 
distributions to a sample of 41 EROs, imaged by HST in the optical or near-infrared 
wide-band filters, in order to identify a set of high-z elliptical candidates on a morpho- 
logical basis. The main results of this work can be summarized as follows. 

a) We have determined the fraction of irregular objects among the EROs of our 
sample, amounting to about 15%. These galaxies are the favourite candidates to be the 
hosts of dusty starbursts. 

b) Considering the whole sample, the galaxies characterized by a brightness distri- 
bution typical of local ellipticals are at least 50%. This estimate is based on a simple 
distinction between irregular/exponentials and dc Vaucouleurs profiles; a more accurate 
analysis, allowing for different values of the shape index n, yields a higher fraction of 
ellipticals (70-80%). 

c) We find that a rather small fraction of EROs (15 ~ 30 %) is made up of compact 
objects whose brightness distributions do not resemble the ones of local elliptical galaxies, 
but are better described by a pure exponential law. 

d) Our data suggest that irregular EROs are found predominantly in the field, and 
that they are characterized - on average - by the reddest colors. 

e) We have determined the rest-frame Kormendy relation for a subsample of 6 cluster 
ellipticals, at redshift ~ 1.3. The relation turns out to be brighter than the local one, at 
fixed size, by 1.4 mag in the B band. 

In the near future, we plan to extend this work to larger, complete samples of high 
redshift galaxies, in particular devoting our efforts to obtain a reliable picture of the 
whole population of high-redshift ellipticals. 
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